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Abstract 

Calculations of time-resolved fluorescence anisotropy from a pair of chro- 
mophores coupled by an excitation transfer interaction are presented. For the 
purpose of investigating the effects of nuclear motion on the energy transfer 



Qj ' and anisotropy, an illustrative model is developed that provides each chro- 

J-^ ' mophore with a single intramolecular vibrational mode. Account is taken of 

(^ ■ non-instantaneous excitation and time- and frequency-resolved detection. Ef- 

C^ ' fects of excitation pulse duration, detection window duration and frequency 

resolution, and excitation transfer coupling strength on the time-resolved 
anisotropy are examined in detail. Effects of vibrational relaxation and 
dephasing are also examined using a simplified Redfield description of the 
effects of coupling to a thermal bath. 
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1. INTRODUCTION 

The process of electronic excitation transfer among chromophores has been widely stud- 
ied. [U In photosynthetic antenna systems, efficient excitation transfer among pigment 
molecules is responsible for greatly enhancing the availability of energy for photosynthe- 
sis. Studies of electronic excitation transfer among impurities in molecular crystals have 
examined the effects of the coupling between electronic and vibrational excitations. ||^ With 
the availability of ultrafast lasers it is now possible to study excitation transfer processes 
occurring on picosecond and sub-picosecond timescales, enabling researchers to directly ob- 
serve some of the fastest excitation transfer processes taking place in photosynthetic antenna 
systems and elsewhere. Generally speaking, the efficiency of excitation transfer depends on 
the proximity of chromophores, their spectral properties, their relative orientations, and 
their interactions with the surrounding medium. When the chromophores coupled by exci- 
tation transfer differ in the orientation of their transition dipole moments, fluorescence (or 
stimulated emission or ground state bleaching) can occur with polarization and frequency 
different from those of the exciting fleld. For this reason, measurement of some variety of 
time-resolved anisotropy can often provide information on electronic energy transport. Two 
experimental methods employing linearly polarized light, pump-probe and fluorescence up- 
conversion, have been used recently to obtain time-resolved anisotropies. The time-resolved 
fluorescence upconversion signal is obtained by combining the fluorescence induced by a laser 
pulse with a delayed gate pulse in a non-linear crystal and collecting the sum-frequency signal 
as a function of the delay time between excitation and gate pulses. By setting the polariza- 
tion of the electric fleld of the excitation pulse parallel or perpendicular to that of the gate 
pulse, parallel (5'||(t)) or perpendicular {S±(t)) components of the transient fluorescence are 
measured. The anisotropy is then calculated according to the expression 

In pump-probe anisotropy measurements, the signal is the relative transmittance of a delayed 
probe pulse. Here the pump pulse is polarized parallel or perpendicular to a probe pulse. 



which passes through the sample, rather than to a gate pulse, which does not. 

Using fluorescence upconversion, Xie et al. ^ recently measured time-resolved anisotropy 
from Allophycocyanin (APC) and C-Phycocyanin (C-PC), the pigment-protein complexes 
isolated from antenna systems of photosynthetic cyanobacteria. The structure of the trimeric 
form of C-PC has been investigated with x-ray crystallography, [^ which revealed that 
chromophores on adjacent monomers are located within ~2 nm of one another, whereas 
inter-chromophore distances within individual monomers are ~5 nm. Xie et al. observed 
sub-picosecond decays in the fluorescence anisotropy in APC and C-PC. Earlier experiments 
by Beck and Sauer employing pump-probe techniques showed that anisotropy in APC 
monomers fails to exhibit sub-picosecond time dependence. Xie and co-workers attributed 
the ultrafast decay in trimeric APC and C-PC to energy transfer between chromophores 
on adjacent monomers. Recent femtosecond pump-probe experiments by Gillbro et al. 
have also observed a sub-picosecond anisotropy decay in C-PC. Rapid anisotropy decays 
have also been observed by Kim et al. |p and Zhu et al. Q in bichromophoric compounds 
9,9'-bifluorene and 2,2'-binaphthyl in solution using fluorescence upconversion |Q and pump- 
probe techniques. [P,PH 

Xie and co-workers [Q made an additional interesting observation that the asymptotic 
value of the anisotropy following the ultrafast decay was sensitive to the excitation and 
detection conditions used. In both C-PC and APC, excitation at the approximate peak of 
the APC and C-PC absorption spectra resulted in signiflcantly lower long time values of 
the anisotropy than excitation closer to the red edge of the absorption spectrum, suggesting 
that a larger fraction of excited chromophores transfer their excitation prior to fluorescence 
in the former case. 

Time-resolved anisotropy can also be a measure of coherence among exciton states, as 
shown in a recent experimental study by Galli et al. |Tl|, where an initial anisotropy of 0.7 
was found in a pump-probe experiment on magnesium tetraphenylporphyrin in solution. 
This flnding supported earlier theoretical work by Wynne and Hochstrasser |T^ in which an 



initial anisotropy of 0.7 was predicted in a molecule having two degenerate electronic transi- 



tions originating from a single ground state, with perpendicular transition dipole moments. 



Wynne and Hochstrasser ||T2[ and Knox and Giilen |T3| have also examined the case of two 
chromophores with a fixed relative angle between their transition dipole moments, showing 
that such a system can exhibit an anisotropy larger than 0.4, the initial anisotropy observed 
in a collection of randomly oriented single chromophores. 

Electronic energy transfer received a thorough theoretical treatment by Forster [|l4l, who 
examined the effects of coupling strength, site energy differences, and details of potential 
energy surfaces and nuclear vibrations on excitation transfer between a pair of chromophores. 
Later work by Rackovsky and Silbey focused on the effects of exciton-phonon interactions 
on excitation transfer among chromophores embedded in a solid. A theoretical study 
by Rahman, Knox and Kenkre [|15[ investigated the fluorescence depolarization resulting 
from electronic energy transfer, pointing out that electronically off-diagonal density matrix 
elements need to be included in order to correctly calculate the polarized fluorescence signal. 

Previous theoretical treatments of time-resolved polarized emission from chromophore 
pairs P, p!B|jr5| (and single molecules with degenerate electronic transitions |]I2|) have ne- 
glected the nuclear degrees of freedom, treating each chromophore as an electronic 2-level 
system. Because femtosecond pulses can launch vibrational wave packets in both ground 
and excited electronic states |16|, the effects of vibrational motion should be included in 



the calculation of time-resolved signals. In terms of the eigenstates of the uncoupled chro- 
mophore pair expressed as products of eigenstates of individual chromophores, we can say 
that Franck-Condon factors and energy level spacings between vibronic states modify the 



effect of a given electronic energy transfer interaction. [Tj] In addition, different sets of states 
can be selected depending on excitation and detection center frequencies and bandwidths. 
In connection to the general issue of simultaneous vibrational coherence and electronic 



energy transfer, it should be pointed out that vibrational quantum beats ^7j and exciton 
dynamics [|I8| have been observed in ultrafast pump-probe anisotropy measurements on 
the light-harvesting complexes of Rhodobacter sphaeroides at 4 K. While an interesting 
argument was advanced |jl8| for the observed spectral and polarization dynamics, no attempt 



was reported to incorporate both nuclear and electronic dynamics in a single model. The 
light-harvesting complex LHl of Rhodobacter sphaeroides has also been the subject of room 
temperature femtosecond fluorescence anisotropy measurements. |T9[ 



Here we investigate a model chromophore pair interacting with a simple detection appara- 
tus in order to explore the roles that molecular parameters and the characteristics of the exci- 
tation and detection processes play on the calculated time-resolved fluorescence anisotropy. 
In order to study the effects of vibrational motion on the time-resolved anisotropy, each 
chromophore is considered to have a single intramolecular vibrational mode. This model is 
presented in Section 2. 

An analytic treatment of a chromophore pair having a minimal number of states, pre- 
sented in Section 3, highlights an important way that chromophore vibrational structure 



can influence the initial anisotropy. Recapitulating the argument of Knox and Giilen |lT3|jr5| 
it is first illustrated how interference in emission from singly excited states with different 
chromophores excited can lead to an initial anisotropy greater than 0.4. The excited states 
must, of course, emit to the same vibronic level of the electronic ground state in order for 
interference to occur. It is pointed out that some vibronic levels are accessible from singly 
excited states having only one specific chromophore excited, and that emission into vibronic 
levels that reveal the identity of the excited chromophore tend to suppress high initial values 
of the anisotropy. 

Having included an intramolecular vibrational mode on each chromophore, we are able 
to examine the roles of simultaneous coherent vibrational and electronic motion induced by 
the ultrashort excitation laser pulse on the anisotropy. Section 4 details these numerical 
calculations. We first choose the parameters of the chromophore pair and the excitation 
and detection processes to correspond as closely as possible to the known parameters in 
experiments by Xie et al. ^ and compare our findings with their experimental results. 
Then we investigate the effects on the anisotropy of varying the molecular and detector 
parameters, such as the strength of the excitation transfer coupling, the difference between 
the zero-zero electronic transition frequencies of two chromophores, and the duration of the 



detection window. In addition to determining time-resolved anisotropies, we also follow 
the time evolution of the excitation probabilities on the two chromophores. Comparing the 
difference in site excitation probability as a function of time to the time-resolved anisotropy 
gives some indication of how the net population transfer affects the observed anisotropy 
signal. 

In Section 5 we examine the effects of vibronic relaxation and dephasing resulting from 
the coupling of the intra-chromophore vibrations to a thermal bath. A simplified Redfield 
theory is used to include these effects. The bath is seen to have a significant influence on the 
time-resolved anisotropy, causing the anisotropy oscillations due to coherent vibrational and 
excitonic motion to diminish in amplitude. It is found, however, that coherent oscillations 
can survive the introduction of fairly rapid vibrational relaxation. Section 6 summarizes our 
conclusions. 

2. MODEL SYSTEM 

We start this Section by describing the model chromophore pair. Each chromophore 
has a ground and an excited electronic state and a single vibrational degree of freedom. 
The minima of the potential energy curves for the ground and excited electronic states 
are spatially displaced from each other in the vibrational coordinate. The Hamiltonian for 
chromophore 1 has the form 

Hi = \gl)Hgi{gl\ + \el){Hei + ei){el\. (2) 

where \gl) and |el) denote the ground and excited electronic states, Hgi and H^i denote the 
corresponding nuclear Hamiltonians, and Si is the zero-zero electronic transition frequency. 
pO| We take h = 1 throughout this paper and thereby express energy in angular frequency 
units. The Hamiltonian for chromophore 2 has a form analogous to Eq. (0). 

The chromophore pair has four electronic states: \gg) (= \glg2)), \eg), \ge), and |ee). We 
omit the site index on the electronic state except when not doing so would create ambiguity. 
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In \gg)., both chromophores are in their respective ground electronic states, \eg) and \ge) 
are the singly excited states with chromophore 1 and chromophore 2 excited, respectively, 
and |ee) is the doubly excited state. The doubly excited state does not participate in the 
electronic excitation transfer and is omitted in what follows, but could come into play in 
pump-probe measurements, where it can be used as the final (measuring) state. [|T^ 



The two singly excited states are coupled by an excitation transfer interaction, which in 
the dipole-dipole approximation would have the form [|14| 



y = ^[P'l- ^^2- ^{lli-'R){ll2-m = U{\eg){ge\ + h.c.), (3) 



where 



ill = lUiGi (|el)(5fl| + h.c.) , 

fl2 = ^2e2i\e2){g2\ + h.c.), (4) 

are the electronic dipole moment operators for chromophores 1 and 2, with fii and fi2 
denoting the magnitudes of the electronic transition dipole moments and ei and 62 denoting 
the unit vectors in the direction of the electronic transition dipole moments of chromophores 
1 and 2, respectively. In the last expression of Eq. @, only the terms that nearly conserve 
electronic energy have been retained, and U is given by 

1 3 

^ = ^3 '"I '"s [(ei ■ es) - -2 (ei • Ro)(e2 ■ Ro)] . (5) 

Combining Eqs. (2-5), the Hamiltonian for the chromophore pair can be written as 

H,p = \gg)[Hgi + Hg2\{gg\ + \eg)[Hei + Hg2 + ei]{eg\ + \ge)[Hgi + H,2 + e2]{ge\ + V . (6) 

The notation we shall use to label the vibrational states of the chromophore pair in the site 
representation refers to each mode's vibrational quantum number with a subscript designat- 
ing the electronic state. For example, \ge)\mgne) is the state with chromophore 2 excited, 
m quanta of vibration in the ground state of chromophore 1 and n quanta of vibration in 
the excited state of chromophore 2. 



The interaction of the chromophore pair with the laser pulse is governed in the dipole 
approximation by 

Ve{t) = -f2cp-E{t) , (7) 

where fl cp = fi i + fi 2, E(t) is the electric field of the laser pulse, and /ii and fl2 are the 
electronic dipole moment operators given by Eq. (|^). We assume E(t) to be a Gaussian 
laser pulse of the form 

E(t) = Eo gl exp (-tV2r^) cos Vtt , (8) 

with amphtude Eq, width r^, related to the FWHM intensity by 

FWHM 

^' = ^TteT- '"' 

center frequency VL and polarization gl. 

The detector consists of a collection of two-level systems with a transition frequency 
distribution taken to be Gaussian. The Hamiltonian for the z*'* two-level detector, which is 
meant to play a role analogous to that of a specific mode of the quantized radiation field, is 
given by 

Hl=\d)El{a'\ + \U)Eim, (10) 

where |a*) and \V) are the ground and excited states of the i*^ detector two-level system, 
respectively, and e^^ = El — E^ is the corresponding transition frequency. The electronic 
dipole operator for the i*^ two-level system has the form 

t2^, = ^^,e,{\a')m + \b^){a^\), (11) 

where /i^ is the magnitude of the transition dipole moment and e^^ is the unit vector specifying 
the detector orientation. We omit the index i on fi^ and e^^, taking all detector two-level 
systems to have the same transition moment regardless of the transition frequency. The 
interaction Hamiltonian of the chromophore pair with the i*^ detector two-level system is 
taken to be 
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\/i = -r^Aiep-/ij,, (12) 

where ?] is a constant that will cancel out in the anisotropy. We can make the rotating wave 
approximation in the chromophore pair-detector interaction, which amounts to only keeping 
those terms that nearly conserve energy. In that approximation, Eq. (O) can be written as 



Vd = -Vf^i f^d (ei ■ Bd) (|6*)|^^)(e^|(a^| + h.c.) 

-7] fi2 /id (62 ■ e^) {\b') \gg) {ge\ {a'\ + h.c). (13) 

Time-dependent perturbation theory is used to determine the state of the system (chro- 
mophore pair + detector) following an interaction with a laser pulse and an interaction 
between the chromophore pair and the detector. Starting in an initial state \gg)\mgng)\a^), 
in which the chromophore pair and the detector two-level system are in their respective 
ground electronic states, the probability amplitude for finding the chromophore pair in the 
vibrational state \m' n') of the ground electronic state and the detector excited to |6*) is 
given in the interaction picture by 

mm'gn'g\{ggm)) = - r dr T dr' 

J — oo J—oo 

xm{m'gn'g\{gg\K{T)UT')\gg)\mgng)\a^) , (14) 

where Ve{T) and V^{t) are the previously defined interaction Hamiltonians, given in Eq. (|^ 
and Eq. (|T3p, respectively, transformed to the interaction picture, 

V,{t) = e'^oWe{t)e-'^ot ^ (15) 

and 

|§(t)) = e'^o^|^(t)), (17) 

with 



Ho = Hep + H^d- (18) 

Notice that the reverse ordering, in which V]_{t') acts before Ve{T) is excluded from Eq. 
(14) under the rotating wave approximation of Eq. ([T3|). The probabihty of finding the 
chromophore pair in the electronic ground state level \Tn' n' ) and the detector in |6*) is then 



given by the absolute value squared of Eq. (|T3): 



nt) = m{my^\{gg\^m'- (19) 

In order to calculate Pn(t) and P]_(t), we set e^^ parallel and perpendicular, respectively, 
to bl, and average over all possible orientations of the chromophore pair in the laboratory 
frame. For given initial and final states of the chromophore pair, we must sum P*(t) over 
the distribution of the detector frequencies. Taking that distribution to be Gaussian and 
converting the sum over detector frequencies to an integral, the total excited state population 
of the detector for given initial and final chromophore pair states becomes 

/oo 
dedD{ed-Sd)Pit,ed), (20) 

-oo 

where 

D{ed - Sd) = -rV exp {-{sd - ^d) V^'} , (21) 

where Ed is the center of of the detector frequency distribution, A is its width, and -P*(t) has 
been rewritten as a function of Ed- 

The total detector population change is obtained by summing over all final states of the 
chromophore pair and a thermal distribution of initial states. The time-resolved anisotropy 
is calculated according to Eq. (P, where the observed emission rate is obtained from the 
time-dependent population of the detector using the expression 

^\\{±){t) - ^^ > (22) 

in which At is the response time of the detector. By calculating the average rate of change 
in detector excited state population, as in Eq. (P^, we crudely incorporate the finite time 
resolution of detection. [^ 
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3. SIMPLE ANALYTIC MODEL 

In Section 2 we developed a model Hamiltonian with which to calculate time-resolved 
anisotropy for a pair of chromophores having multiple vibrational levels and taking into 
account time and frequency resolution effects encountered in experimental measurements. 
Before proceeding to examine these effects, we shall first investigate the properties of a far 
simpler system, ignoring coherent vibrational motion by keeping only the lowest vibrational 
level in each of the the singly excited electronic states, \eg)\OeOg) and |5'e)|0g0e), and the three 
lowest vibrational levels of the ground electronic state, \gg)\OgOg), \gg)\lgOg), and \gg)\Oglg). 
In their recent work, Xie [^, Wynne and Hochstrasser |jl2|, and Knox and Giilen |ll 



considered similar model systems, determining that if the initial state (prior to excitation) 
and final state (following spontaneous or stimulated emission) is \gg)\OgOg), then the initial 
anisotropy can take on values greater than 0.4. The theoretical maximum value of 0.7 is 
exhibited by pairs of chromophores with perpendicular transition dipole moments of equal 
magnitude. Both Wynne and Hochstrasser |T^ and Knox and Giilen [|13[ have included 



a phenomenological treatment of relaxation between the excited states of the system. We 
ignore relaxation in this section, but the analytical results for the time-resolved anisotropy 
of the simple model considered here will illustrate the dependence of the anisotropy on the 
final state(s) of the system. 

The Hamiltonian of the simplified chromophore pair becomes 

Hep = \gg){gg\{Egg\OgOg){OgOg\ + {Egg + ^ ) | 1 ,0, ) ( 1, 0, | + {Egg + U) \0 gl g) {O gl g\} 

+ \eg){eg\Eeg\OeOg){OeOg\ + \ge){ge\Ege\OgOe){OgOe\ + U{\eg){ge\ \0e0g){0g0e\ + h.c.) , (23) 
Diagonalizing the excited state portion of H^p yields the exciton states of the system, given 

by a 

n n 

1+) = sin -15-6)10306) + cos-|e5()|0e0g) , 

n n 

I-) = cos-|5e)|0,0e) -sin-|ef7)|OeO,) , (24) 
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where tan^ is 



ta„»=?^M^, (25) 

and the eigenenergies are 

E± = ^{E,g + Eg,) ± \^{E,g - E g.f + 4f/2 ( ^ 0, | 0, 0^ ) ^ . (26) 

The form of the dipole moment operator in the exciton representation becomes 

/iep = M+{|^^)|OA)(+l + h.c.} + M-{|^^)|OA)(-| + h.c.} 

n n 

+ ^ii{lgOg\OeOg) {cos -\gg)\lgOg){+\ - sin -\gg)\lgOg){-\ + h.c.} 

n n 

+fi2{0glg\0g0,){sm-\gg)\0glg){+\+cos-\gg)\0glg){-\ + h.c.} , (27) 

with n + and /J, _ given by 

n n 

^+ = H2{0g0g\0g0e)sm- +/^i(0g0g|0e0<,)cOS-, 

n n 

fl =^l2{0g0g\0g0e)cOS--f^l{QgQg\Qe0g)sm-. (28) 

In addition to simphfying the chromophore pair model, we also simplify the interaction 
with the laser pulse and the detector by assuming that the excitation pulse duration and 
the detection window are both much shorter than the timescales of dynamics in the excited 
electronic state. Under these conditions, the rate of emission, S(t), can be obtained by 
differentiating both sides of Eq. (^) with respect to time, 

/oo 
dedD{ed~ed)Sit,ed), (29) 

-oo 

where S{t,ed) is the time derivative of P(t,ed) given by Eq. (|I^) and D{ed — Sd) is defined 
in Eq. (^11) . The probability amplitude for finding the system in \gg)\QgQg)\b), starting in 
the initial state \gg)\QgQg)\a), is given in the interaction picture by 

~ 1 /■* 

{h\{QgQg\{gg\^{t)) = — -=V7rrLEo?7yU^/id(e+ ■ eL)(e+ ■ e^) J drexpi^Ed - e+g)r 

1 /■* 

— y=VnTLEoT]fL^_fid{e_ ■ Bl) (e_ ■ e^) J dr expired - e^g)T , (30) 
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where e±g = E± — Egg and n±= /i±e±. Calculating S{t, Ed), integrating over the frequency 



profile of the detector according to Eq. (p9D, and taking A :^ I&+-I where e+_ = E^ — E_ 
yields 

Sit) = \'nrlEW^,X^,l{e^ ■ e,)\e^ ■ e,)^ + \1xrlEW^,t^?,{e^ ■ e,)^(e_ ■ e,)^ 

+ TTTlElrf ^\^i^__^l{e+ ■ eL)(e+ ■ ed)(e_ ■ eL)(e_ ■ e^) cos (e+_t) . (31) 



The first two terms in Eq. (|3lD result from excitation and emission via |+) and |— ), respec- 
tively, and the third (oscillatory) term results from interference in emission from |+) and 
|— ). While the internal geometry of the chromophore pair is assumed to be fixed, it can be 
randomly oriented in space. To average Eq. (|31|) over spatial orientations, we integrate e+ 
over 47r of solid angle and, maintaining a fixed angle between e_|_ and e_ , integrate the latter 
over 27r radians. Defining a direction cosine 

cos 7 = e_|_ ■ e_ , (32) 

the orient at ionally averaged signal is given by (see Appendix) 

S^it) = ^^^^^^{^(/it + -"-) + 2/i>^_(icos% + ^sin^T) cos (e+_t)} , (33) 

Snt) = ^^'^f'^n ^(/^t + /^-) + 2/i>^_(^cos^ - Isin^) cos (.+_t)} . (34) 
Using Eqs. (0) and (Q) in Eq. (|l]) gives the time-dependent anisotropy 



^ 2(^t + /^-) + /^i^-(4 + cos^ 7) cos (£+-t) 
^' 5(/it + ;ui) + 10/iV-Cos2 7cos(£+_t) ■ ^^ 

Let us now examine the case of identical chromophores. In that case, and assuming 
6* = 5- in Eq. (^), the expressions for yU+ and yU_ in Eq. ( pH]) simplify to 

H±e± = ^(0g0g|0e0g)(/i2e2 ± yUiGi) . (36) 



Because /ii = /i2, the two transition dipole moments in Eq. (|36D are perpendicular. Substi- 
tuting 7 = 7- into Eq. ( pS]) yields 
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2(;U^ + /ii) + 3 yu^yu?. COS (e+_t) 
5(/i+ + fit 



^ ^KH'+^H'-) -T o/x^AJ._L-uoi^c+_i.; 



The expression (R^ for the time-resolved anisotropy can be restated in terms of the angle 



between the directions of the transition dipole moments of the two chromophores, ei and 
62- Calling that angle 0, Eq. ( ^7\) becomes (see Appendix) 

R(t) = 0.4 + 0,3 ^ -^QS0 cos (£+_t) . (38) 

1 + cos"' (J) 

The denominator of Eq. ( P7| ) is proportional to the total intensity and is time- 
independent. The time-independent overall intensity results from the fact that the chro- 
mophore pairs are identical and randomly distributed. When transitions to both exciton 
states are optically bright, the oscillatory term initially enhances the parallel signal while 
diminishing the perpendicular signal, resulting in an initial anisotropy larger than 0.4. Al- 
ternatively, if one of the exciton states is optically dark, the interference term vanishes and 
the anisotropy is 0.4 and time- independent. The maximum anisotropy can be attained when 
the interference term is maximized, which occurs when ei and e2 are perpendicular. 

Let us now compare the anisotropy in Eq. (|38| ) with the anisotropy that would be 

obtained if the only admitted final states of the chromophore pair were \gg)\lgOg) and 

\gg)\Oglg), the two electronic ground states with one of the chromophores in the first vi- 

brationally excited state. The probability amplitude of finding the system in \gg)\lgOg)\b), 

starting in the initial state \gg)\OgOg)\a) with an excitation pulse short on the timescale of 

excited state dynamics, is 

~ 1 f) 

{b\{^gOg\{gg\'^{t)) = — ■y=^/^TLEor|f^dCos-{OeOg\lgOg)fi+fil{e+■eL){el■ed) 

rt 
X 



Jo 



1 f) 

+ ^v^rL£'o?7/irfSin-(OeOg|lgOg)/i__/ii(e_ ■ ei)(ei ■ e^) 



X 



/ dTexpi{ed — £-g + uj)t . (39) 



Since the Franck-Condon factor between states |OgOe) and |lgOg) vanishes, there is no e2 
component present in emission from |+) and |— ) into \gg)\lgOg). The rate of emission, 
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obtained in the same manner as Eq. (|3TD , is 



S{t) = ^TTT^E^^V+Zii/idCos^ ^(0,0,|l,0g)2 (e+ ■ ei)2(ei ■ e,y 



.2 7-i2„2 2,2 ■ ^ ^ 



+TTTLEQr] /i+/i_/ii/i^ sin- COS -(OeOg I IgOg) (e+-eL)(e_ ■ eL)(ei ■ 6^)2 cos (e+.t) , (40) 

where the overbar signifies orientational averaging. The rate of emission into \gg)\Oglg) 
has an identical form to Eq. (PPJ). The time- resolved anisotropy for the case of identical 
chromophores now takes the form 

r,.,. 2(/it + ^i) + /iV,(6 COS (e+_t) - 2) 

in terms of the transition moments of the exciton states or, alternatively, 

R{t) = 0.1 + 0.3 cos^ + 0.3(1 - cos^ 0) cos (£+_t) , (42) 

where is the angle between ei and e2. The time zero anisotropy predicted by Equations 
( JT) and (J2) is 0.4, independent of the angle between ei and 62. 



The quahtative difference between anisotropies described by Eqs. (|38D and (^) at 
t = can be understood by considering Eqs. (|3T|) and (|40|) in the short time limit where 
cos{e +^6t) = 1. The emission rate in Eq. (^Tj) at 6t becomes 

Sm ^ KrlEyfitfiliOgOglOeOgf (ei • eL^ie, ■ e^)^ 



+ -7TTlE^,r]^f,tf^l{0g0g\0g0,)^ {b^ ' B L^io^ ■ erf)2 



+7rr2EoV/^?/^2/^d(OgOg|Oe03)(OgOj030e) (ei ■ ei)(e2 ■ ei)(ei ■ ed)(e2 ■ e^) , (43) 

where we have used Eq. (^) to express the rate in terms of the site transition dipole 



moments. Carrying out the orientational average (see Eq. (88) of the Appendix) and 
calculating the the anisotropy yields 

^ 2fif{0g0g\0e0g)' + 2/i^ (0,0g 1 0,0^ ) ^ + /^^/^^(cos^ <P + 3 ) ( OgOg 1 0^0,) (0,0g|0,0e) 
^ ' 5/if(O,O,|OeO,)2 + 5/x|(O5O,|O3Oe)2 + lO/xf/xicos2 0(O,O,|OeO,)(O,O,|O,Oe) ' ^ ' 
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where is the angle between ei and 62. The first two terms in the numerator and the 



denominator of Eq. ( ^^ arise from contributions involving each chromophore separately, 
while the third term results from contributions from both chromophores. This crossterm 
depends on the angle between the dipole moments of the chromophores, and the short-time 
anisotropy has its maximum when = S. On the other hand, the rate of emission into 
\gg)\lgOg), given in Eq. (^OD , becomes 

S{St) = ^7rTlEyfj,tfj,l{0g0g\0e0g){0e0g\lg0g) (gi ■ eL)2(ei ■ e,)2 , (45) 

where Eq. (^) was again used to express the rate in terms of transition dipole moments 
of the individual monomers. Carrying out the orientational average and calculating the 



anisotropy yields R{6t) = 0.4. There is no crossterm in Eq. (^Sf ) analogous to the one in 



Eq. ( ^31) because transitions from \eg)\OeOg) to \gg)\Oglg) and from |5fe)|0c,0e) to \gg)\lgOg) 
are Franck-Condon forbidden. There is only one route to each final state in this case, and 
quantum mechanical interference in emission from the two sites does not occur. Therefore, 
the t = anisotropy given by Eq. ( |5BD is sensitive to the presence of both chromophores, 
while the anisotropy given in Eq. (^2|) is not. It is worth emphasizing that the high initial 
value of the anisotropy predicted by Eq. ( PD reflects the fixed relative geometry of the 
chromophore pair, rather than the presence of energy transfer, per se. The lack of dependence 
of S{6t) in Eqs. (^) and (^51) on the exciton splitting is in keeping with that fact. 

Although it might appear in Eq. (|38D as though R{t) would exhibit oscillations whenever 
7^ 0, this is not always the case. For example, in the case when = S and either ei or 62 
lies along the inter-chromophore axis, the energy splitting between the exciton states will 
be zero, assuming dipolar coupling of the form given in Eq. (j^), and the anisotropy will be 
time-independent. Conversely, excitation transfer is not always required for the anisotropy 
to exhibit oscillations. The anisotropy given by Eq. (|35|), which applies to chromophores 
that are not necessarily identical, will exhibit oscillations even in the case of zero excitation 
transfer coupling, as long as the individual chromophores have different electronic transition 
frequencies. 
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From the calculation presented in this Section we see that time-resolved anisotropy in 
a multilevel system will contain both kinds of contributions described by Eqs. ( pSf ) and 
(^), and that proper treatment of time- resolved fluorescence anisotropy must include the 
multilevel structure of both the excited and ground electronic states of the chromophores. 
The relative weights of the contributions leading to the anisotropics given in Eqs. ( P5| ) and 
(^) will be determined by the spectral characteristics of the excitation and detection pro- 
cesses. Generally speaking, contributions giving rise to the standard (0.4) value of the initial 
anisotropy are expected to outweigh the contributions of the kind described by Eq. (^8|), 
especially when higher vibrational states of the ground electronic state are selected as the 
final states by the detection process, which is often the case with fluorescence upconversion. 



4. EFFECTS OF VIBRATIONAL MOTION AND EXCITATION/DETECTION 

CONDITIONS 

The analysis of Section 3, despite using a simplified model chromophore pair and neglect- 
ing effects of less-than-ideal excitation and detection conditions, serves to demonstrate the 
importance of properly including the vibrational structure of initial and final states when 
calculating the fluorescence anisotropy. In this Section, we will utilize the model developed 
in Section 2 more fully, including the full vibrational level structure of both chromophores 
and non-zero duration for the excitation pulse and the detection window. 

The model in Section 2 was developed in the site representation in order to give a better 
physical picture of different features of the model. The calculations in this Section, like those 
in Section 3, are more conveniently carried out in the energy eigenstate representation. The 
chromophore pair Hamiltonian in the eigenstate representation can be written as 

H,, = Y.\a)E^{a\+Y.\^)E,{^\ (46) 

a 7 

where, for simplicity of notation, we denote all eigenstates in the ground electronic manifold 
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of the chromophore pair by a and all exciton states by 7, with E^ and E^ as the correspond- 
ing angular frequencies. The electronic dipole operator can be expressed as 

/icp = ^/^a7ea7(|a)(7l + l7)(a|) , (47) 

a, 7 

where iia-y^a-y is the transition dipole moment between states |a) and I7). 

Starting in the initial state |a)|a), the probability amplitude for finding the system in 
a final state | «')!&) following an interaction with a laser pulse and the detector is given by 
(see Eq. (0)) 

a'\^{t)) = ^^7— ^Xl/^"7/^"'7(^a7 ■ ei)(ea'7 " Grf) / dr j (Iti 

Z ^, J— 00 J— 00 



X 



^ J— 00 J— 00 

exp {i{€d - €ya')r} exp {i{6-ya - fi)ri} exp {-T^/2rl} . (48) 



Calculating the excited state population in \b) according to Eq. (|T9D and integrating over 
the frequency profile of the detector according to Eq. (|20D yields 

n2^2,,2 



7:7' 

X / dr dri dT2 / (iracxp {-(r - r2)^A^/4} 

J— 00 J— 00 J —00 J~oo 

X exp {-i{ed - eya')^} exp {-i(£:y„ - fi)ri} exp {-t^/2tI} 



X 



exp {i{ed - £700^2} exp {i{e^a - ^)t"3} exp {-t^/2ti} (49) 



As discussed in Section 2, the total excited state population is calculated by summing Eq. 
(^) over all final states and thermally weighted initial states of the chromophore pair. 

The average, over chromophore pair orientations, of the products of four direction cosines 
appearing in Eq. ( ^91 ) must now be carried out. It is convenient to express the transition 
moments e^^yfia-y in terms of ei/ii and 62/12, the transition dipole moments of the individual 
chromophores, 

ea7/ia7 = {a\gg) {eg\-f)eifii + {a\gg){ge\-f)e2iJ,2 ■ (50) 



Using Eq. ( |5DD enables us to express the product of four direction cosines in Eq. 

in terms of direction cosines involving ei and 62- There are four different orientational 



averages to perform. These averages are given in the Appendix. Eq. (^9]) can now be 
used to compute the population in the excited state of the detector, and subsequently the 
time-resolved fluorescence anisotropy using Eq. (p2D. 

In the numerical calculation presented in this Section, we have chosen the parameters 
for the chromophore pair to roughly correspond to the known parameters in the closely- 
coupled chromophore pairs consisting of a-8A and (3-8A tetrapyrroles on adjacent monomers 
in C-PC trimers. @] The value of the excitation transfer coupling matrix element is set to 
-50 cm^^ (1 cm~^ corresponds to an angular frequency of 1.885 x lO^^s^^), corresponding to 
Sauer and Scheer's calculation of the excitonic interaction energy. p4| The angle between the 



transition moments of the chromophores is set to 65°, corresponding to the value estimated 
from the x-ray crystallographic structure. The zero-zero electronic transition frequency 
of chromophore 2 is taken to be 150 cm"^ higher than that of chromophore 1. This offset 
corresponds to the difference in absorption maxima of the a-84 and (3-84 chromophores. 
p^ , |25[] The potential energy curves for Hgi, Hg2, H^i, and He2 are chosen to be harmonic and 
to have the same frequency of 100 cni"^. This choice of vibrational frequency is arbitrary, but 
is typical of a low frequency molecular vibration and is in the interesting range comparable to 
the frequencies of pure exciton dynamics. In the absence of experimental data, we assume a 
small displacement between ground and excited potential energy curves on each chromophore 
of 1.5 times the rms width of the ground state coordinate distribution, -^ for a mass- 
weighted oscillator coordinate. The excitation laser pulse is taken to be 70 fs Full Width at 
Half Maximum (FWHM) in intensity {t^ = 42 fs), in accordance with experiment. |^ The 
frequency resolution of the detector is set to 210 cm~^ FWHM in intensity (A = 126 cm~^), 
which corresponds to the frequency width of the detection window used in the experiments, 
and the width of the time detection window is taken to be 70 fs. The radiative lifetimes 
of the excited states are much longer than our timescale of interest of a few picoseconds, and 
non-radiative excited state population decay mechanisms, such as intersystem crossing, are 
ignored. The initial state of the chromophore pair is \gg)\OgOg) for all calculations in this 
section. See Section 5 for calculations that include additional thermally populated initial 
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states. 

Figure 1 shows time-dependent anisotropies calculated for the model system using the 
parameters listed above. Two different sets of excitation and detection conditions are illus- 
trated. The upper (thick) curve in Figure 1 was calculated with Q, the center frequency of the 
laser pulse, resonant with ei, the zero-zero electronic transition frequency of chromophore 1. 
The lower (thin) curve was calculated with Q at 100 cm~^ above £2, the zero- zero electronic 
transition frequency of the chromophore 2. In each case, the average detector frequency, Ed, 
was centered 250 cm~^ below Q. We included 30 states in the singly excited manifold, allow- 
ing up to 4 quanta of combined vibrational excitation, and 28 states in the ground manifold, 
allowing up to 6 quanta of combined vibrational excitation. Inclusion of additional states in 
the ground or the excited manifold had no effect on the calculated anisotropy. 

The difference in excitation and detection conditions in the two cases gives rise to qualita- 
tive differences in the anisotropy. Anisotropies in Figure 1 exhibit complicated behavior that 
contains effects of both exciton dynamics (including excitation transfer) and nuclear dynam- 
ics. It can be seen from the near constancy of the anisotropy that little change occurs in the 
emission polarization direction when the excitation pulse only weakly populates the excited 
states of the higher-energy chromophore (thick curve). The absence of exciton dynamics is 
reflected in the weak time-dependence and relatively high value of the anisotropy. On the 
other hand, anisotropy resulting from excitation at 100 cm~^ above the zero- zero electronic 
transition frequency of the higher-energy chromophore, shown by the lower (thin) curve in 
Figure 1, resonantly populates a range of states having amplitude on both chromophores, 
leading to a larger extent of excitation transfer, as evidenced by the large oscillations in the 
anisotropy. The dependence of the anisotropy on the excitation frequency in our calculations 
shows a similar trend to the experimental results of Xie et al. [Q , who found that excitation 
at the red edge of the absorption spectrum of the chromophore pairs in APC and C-PC 
trimers resulted in a more weakly time-dependent anisotropy with higher asymptotic values 
than excitation at a higher frequency. 

One feature common to the two plots in Figure 1 is their high initial anisotropy and 
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subsequent rapid decrease during the first 0.1 ps. Despite the difference in site energy 
between the chromophores and the finite bandwidth of the excitation pulse and the detector, 
both chromophores can be excited at short times and the interference in emission necessary 
for anisotropy greater than 0.4 can occur. When the excitation pulse and the detection 
window overlap, both are effectively shortened by the requirement that the former must act 
before the latter. In addition, the nuclear wave packets \ipeg{t)) and {ipgeit)) prepared in the 
Franck-Condon regions of the two excited state potential energy surfaces both resemble \OgOg) 
at short times, favoring emission into the same ground state. Thus, the initial anisotropy is 
dominated by terms of the kind illustrated by Eq. (|38D . 

The effective frequency resolution improves once the pulses cease to overlap and each can 
operate for its full duration. Moreover, since \ipeg{t)) and \ipge{t)) move differently on their 
respective potential energy surfaces, nuclear motion favors emission to non-overlapping final 
states, decreasing the effect of quantum mechanical interference in the anisotropy. It should 
also be mentioned that the site energy difference between the two chromophores contributes 
to a differing rate of phase accumulation, which leads to to an overall cosinusoidal oscillation 
in the anisotropy. All of these factors can contribute to a rapid decrease in the anisotropy, 
even in the absence of energy transfer coupling between the sites. For this reason, the rapid 
"decay" of the anisotropics of Figure 1 cannot automatically be regarded as a manifestation 
of ultrafast energy transfer. 

In the calculations presented in Figure 1, the time detection window At was taken to be 
70 fs, equal to the FWHM duration of the excitation pulse and the time corresponding to the 
FWHM frequency width of the detector. However, experimental time resolution can often 
be lower than the inverse frequency width of the detector. Figure 2 compares the anisotropy 
shown in the lower (thin) curve of Figure 1 with the anisotropy calculated using identical 
parameters except that the time detection window. At, was lengthened to 150 fs. The t = 
anisotropy drops considerably as the time resolution decreases. Therefore, limitations in the 
experimental time resolution can prevent observation of the high initial anisotropy discussed 
above. 
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Figure 3 shows time-dependent anisotropies calculated with the same parameters as 
those used in Figure 1 except that the offset between the zero-zero electronic transition 
frequencies of the two chromophores is decreased to 100 cm~^ (i.e. equal to the vibrational 
level spacing) by decreasing the zero- zero transition frequency e^ of the higher-energy chro- 
mophore, resulting in more efficient excitation transfer because many vibronic levels of the 
two chromophores are now in resonance. Because the excitation frequencies were kept the 
same as those in Figure 1, the fi value for the lower (thin) curve in Figure 3 is now 150 
cm~^ above the zero- zero electronic transition frequency of the higher-energy chromophore. 
The time-resolved anisotropies in Figure 3 show considerable differences from those of Fig- 
ure 1. The lower (thin) curve in Figure 3 exhibits more frequent large oscillations, due to 
more efficient energy transfer between the singly excited states. The upper (thick) curve in 
Figure 3 shows an increase, relative to the corresponding case of Figure 1, in the magni- 
tude of oscillation, consistent with the fact that there is now more spectral overlap of the 
excitation pulse with states having significant amplitude on the higher-energy as well as the 
lower-energy chromophore. 

In Figure 4, time-resolved anisotropies are shown for a case of two identical chromophores, 
with the same excitation and detection conditions as in the previous calculations. The 
anisotropies can be seen to refiect the fact that this is the case with the most efficient 
excitation transfer. The frequency of the large amplitude oscillations is the highest compared 
to Figures 1 and 3 in the case of excitation 250 cm~^ above the zero-zero electronic transition 
frequency of the chromophores (thin curve). Excitation at the zero-zero electronic transition 
frequency of the chromophores (thick curve) shows the largest extent of excitation transfer 
(even leading to a negative anisotropy at t ~ 330 fs) compared to its counterparts in Figures 
1 and 3. 

Figure 5 examines the dependence of the initial it = 0) anisotropy on the strength of 
the excitation transfer coupling while the rest of the parameters are kept constant. As the 
excitation transfer coupling is increased, the initial anisotropy decreases because the initial 
rapid decrease of the anisotropy, calculated before time averaging according to Eq. (p2[) 

22 



(not shown), becomes more rapid with increasing couphng, and keeping the time-resolution 
of the detector constant results in observing a lower average anisotropy during the t = 
detection window. It is evident that for the case where the timescale of exciton dynamics is 
shorter than the timescale of the excitation and detection process, the high initial anisotropy 
will not be observed. 

The time-resolved anisotropy reflects the exciton dynamics of the chromophore pair. It 
would be interesting to determine the extent to which the time-resolved anisotropy mani- 
fests the net transfer of excited state population between the chromophores. [^ The time- 
dependent site population difference is given by (\I'(t)|Pi — P2\'^{t)), where |\I'(t)) is the 
state of the chromophore pair at time t, and Pi an P2 are given by 

Pi= Yj \eg) {eg\ In^mg) {nerugl (51) 

ricirig 

and 

P2= Yj \ge){ge\\ngme){ngme\ . (52) 

ng,me 

At times after the excitation pulse has subsided, the time-dependent orient at ionally averaged 
population difference can be expressed by 



(^(t)|Pi - P2|*(i)) = X!/^"7/^"7'(ea7 ■ eL)iea^' ■ e^) 



7,7 



xF{e,^ - Q)F{ey^ - fi)e^^-v*(7|(p, - p,)\^') , (53) 



where 



F{uj) = V2nTL exp (-r£cj72) . (54) 

Figure 6 shows the calculation of the expectation value of Pi{t) — P2{t), normalized by the 
total excited state population, Pi(t) + P2(t), for the same parameters as in Figure 1. Again, 
the initial state of the chromophore pair is taken to be |a) = \gg)\OgOg). The thick curve 
shown in Figure 6a is the time-dependent site population difference in the case of excitation 
centered at the zero- zero electronic transition frequency of chromophore 1, the lower-energy 
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chromophore, and the thin curve is the corresponding time-dependent anisotropy. In Fig- 
ure 6b, the thick curve is the time- dependent site population difference for the excitation 
frequency centered 100 cm'^ above the zero-zero electronic transition frequency of (the 
higher-energy) chromophore 2, and the thin curve is the again the corresponding time- 
resolved anisotropy. The excitation condition that leads to the site population difference 
curve in Figure 6a places most of the initial excitation on chromophore 1, as evidenced by 
the positive initial value of r}\,l , — t^),1 - The relatively small oscillations of the upper 
(thick) curve indicate that most of the excitation remains on the lower-energy chromophore. 
This is consistent with the anisotropy results shown in the thin curve of Figure 6a, where 
the time-resolved anisotropy obtained with the same excitation conditions shows relatively 
small oscillations and high anisotropy. Excitation centered 100 cm"^ above the zero-zero 
electronic transition frequency of the higher-energy chromophore, resulting in a site pop- 
ulation difference shown in the thick curve of Figure 6b, also exhibits behavior consistent 

P (t) — P (t) 
with the excitation conditions, p /A , p )A is initially negative, indicating that the excita- 
tion conditions place most of the excitation on the higher-energy chromophore; even though 
there is good spectral overlap between the excitation pulse and states on both chromophores, 
the Franck-Condon factors are more favorable for the excitation of the higher-energy chro- 
mophore for the coordinate displacement between ground and excited states that we have 
chosen. The large oscillations about a relatively low value of the site population difference 
indicate efficient net population transfer between the two chromophores. This behavior is 
also reflected in the large oscillations seen in the corresponding time-dependent anisotropy 
shown in the thin curve in Figure 6b. In fact, the flrst several extrema in the population 
difference of Figure 6b coincide with corresponding extrema in the anisotropy. In both Figs. 
6a and 6b, however, the anisotropics show high frequency structure attributable to coherent 
vibrational motion, that is not as clearly evident in the excitation probability difference. 
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5. EFFECTS OF VIBRATIONAL RELAXATION AND DEPHASING 

In the previous sections we developed a model for an isolated chromopliore pair and a 
detector in order to describe the time-resolved fluorescence anisotropy following excitation 
with an ultrashort laser pulse. The time-resolved anisotropies shown in Section 4 exhibited 
a similarity in excitation wavelength dependence to the experimental measurements made 
on APC and C-PC [Q The experimental anisotropy measurements did not, however, show 
evidence of the coherent oscillations in the time-resolved anisotropy that were seen in the 
calculations on isolated chromophore pairs. The surrounding medium undoubtedly plays a 
major role in destroying the coherences between exciton states that give rise to the oscil- 
lations seen in the calculations of Section 4. In this section, we include two aspects of the 
surrounding medium, namely, vibronic energy relaxation and dephasing, and examine their 
effects on the time-resolved anisotropy. 

Our inclusion of vibrational relaxation and dephasing builds upon the approach of Jean, 



Friesner and Fleming. [^ In their work, Jean et al. developed a quantum mechanical theory 
of photo-induced electron transfer (see also the recent work of Pollard and Friesner on the 
application of Redfield theory to multilevel systems [^ and work by Walsh and Coalson 



||28||). Using Redfield relaxation theory, they incorporated vibrational relaxation through a 
coupling of the quantum-mechanical reaction coordinate of an electron transfer system to 
a thermal bath. We adopt a similar approach in our work, introducing a coupling of the 
intra-chromophore vibrations to a thermal bath. 

In Section 5A, we derive the form of the time- and frequency-resolved signal in terms of 
the elements of the reduced chromophore pair density matrix. Then we introduce and de- 
scribe the system-bath coupling operator in Section 5B. Section 5C explains some important 
details of our Redfield calculations, and results of numerical calculations of the time-resolved 
anisotropy including vibrational relaxation and dephasing are presented in Section 5D. 

5 A. Density Matrix Increment 

In deriving a density matrix expression for the time- and frequency-resolved signal, the 
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notation of Section 4 will be used. In order to calculate the excited state population of the 
detector we need the contribution to the density matrix that is second order in the electric 
field of the excitation laser pulse and second order in the interaction with the detector. This 
contribution is given in the interaction picture by 

Jt(-) Jto Jto Jto 

+ V,{T)V,{Ts)poV,{T2)V,{n) + V,{T)Ve{T2)poV,{r,)V,{n) + h.c.} , (55) 

where Vd{t) and Ve{t) are the previously defined interaction Hamiltonians in the interaction 
picture, given by Eqs. (|T^ and (PIB]), and po is the initial density matrix of the entire system, 
which is assumed initially to take the factorizable form 

Po = Pep Pb \a){a\. (56) 

Pep = |q^)(q^| is the initial reduced density matrix of the chromophore pair and pi, is the 
initial density matrix of the thermal bath. Hq now includes the Hamiltonian of the thermal 
bath, Hf,, and the chromophore pair-bath coupling, Hf,cp, and is given by 

Hq = Hep + Hd + Hb + Hbcp ■ (57) 

The precise form of Hbcp will be developed in Section 5B. 

If the time t in Eq. (^) is large enough so that excitation and detection do not overlap, 
and the density matrix increment is defined as Ap^^'^ (t) as 

Ap(^)(t)^p(^)(t+^)-pW(t-^), (58) 

then the detected signal can be calculated with the help of 

Ap(^)(t) = rj dr r dn r dr^ r d^ Vd{T)Ve{T;)p^V,{T2)Vd{n) + h.c. (59) 

Jt-^ Jtn Jtn Jtn 



We arrive at Eq. ( |59D by noticing that when excitation and detection do not overlap, the 
first term on the right-hand side of Eq. (^) does not contribute to Eq. (|5^, and the upper 
limit of integration in the integral over T2 can be replaced by infinity. 
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The integrand in Eq. (|59[) can be simplified as follows. First, we make an assumption 
that the effects of the bath are negligible on the timescale of the duration of the laser pulse. 
That approximation allows us to neglect the system-bath coupling Hamiltonian in the terms 
Vei^Ts) and ^(^2), and Eq. ( |59D becomes 

xF(£^« - fi)F(£y„ - n)Vd{T) |7)(7'| pb \a){a\Vd{n) + h.c. (60) 



where F{e^a — ^) has been previously defined by Eq. (]5^). We next write out explicitly the 
operators Vrf(r) and Vrf(ri) in Eq. ( pDj ) and take the expectation value of Ap*^'')(t) in state 
\h)\a'). Since the Hamiltonian of the detector commutes with the rest of the operators in 
ifo; we replace Hd with the corresponding eigenenergy. Integrating the resulting expression 
over the frequency profile of the detector according to Eq. (pO]) yields 

' dedD{ed-edm{a'\^P^'^\t)\a')\h) = 

D 

7712 2 2 
-, zZ /^a7/^a7'/^a'/3/^a'/3'(6a7 ■ 6i)(6a7' ■ ^i)(^Q'/3 ■ ed)(eQ,/^/ • erf) 



7,7',/9,/3' 



/■-I-— T 

X I J rfr I driF(e^„ - (])F(£y, - (])e'^''("-"^)e-("-"i)'^'/^ 

where D{ed — Sd) denotes the frequency window of the detector, given by Eq. (0), and 
the prime on H'q indicates that Hd has been removed. At this stage we make a similar 
assumption regarding the coupling of the chromophore pair to the bath during detection as 
we made during the excitation pulse, namely, that the bath has a negligible effect on the 
chromophore pair on the timescale of the inverse of the frequency width of the detector, 
A^^. That assumption allows us to make the following approximate replacement in Eq. 

m 

(/5|e-'^0"|7)(7'| Pb e'^0"^|/5') ^ (/3|e"'^o-|^)(y| p, eiH',r^py-i{H,+E,,,){r-r,) _ ^g2) 

Using Eq. (^) in Eq. (^Tj) and taking the trace over the bath results in the final expression 
for the population of the state \a')\h): 
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4 

7,7',/3,/3 



/oo 
rf£,D(£,-5,)(6|(a'|ApW(t)|a')|6)} = 



xTr,{(/3|e-^^o-|7)(7l Pb e^^o-|/5')} + c.c. (63) 



Making approximations that neglect the effects of the bath during the excita- 
tion pulse and the detection time window confines the effects of the bath to 
Trb{(/5|e~'^o'^|7)(7'| pb e^^o'^\P')}. In order to evaluate the latter quantity we must now 
develop a model for the bath and the form of the bath-chromophore pair interaction. 

5B. Bath-Chromophore Pair Interaction 

Before proceeding farther, we make the following definition regarding the system and 
the bath. The degrees of freedom that we choose to treat explicitly are the two intra- 
chromophore modes, referred to as Qi and Q2- In general, we assume that any optically 
active collective coordinates (i.e. those whose equilibrium values are displaced in the excited 
electronic states) would be treated on the same explicit footing as Qi and Q2. Thus the 
bath modes are defined to be those collective coordinates that are not displaced when an 
optical transition takes place. 



Consistently with our assumption of the factorizable form, Eq. (|56D of the total initial 
density matrix in the electronic ground state, we neglect bath-chromophore pair interac- 
tion entirely prior to electronic excitation. We expand the system-bath coupling operator 
VegiQi, Q2, q) for the singly excited state \eg) about the equilibrium positions of Qi and Q2 
in that state, denoted by A^^ and Ag^, respectively, and let q denote the bath coordinates. 
Expanding Veg{Qi, Q2, q) to first order in the displacement of Qi and Q2 yields 

(ft - Arf-M^^. (64) 



The g-dependence of the first term on the right-hand side of Eq. (p4\ ) would, in general, 
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contribute to optical dephasing. Considering g to be a single degree of freedom for purposes 
of illustration, we expand the first term on the right-hand-side of Eq. ( [6^ ) in powers of q 
about its equilibrium position, which is chosen to be zero (and by definition, to be the same 
in the ground and both excited electronic states): 

ye,(Af , A^^ q) - l-,,(Af , A^/, 0) + Y '''^^^, '' ^ + ... (65) 

where the term linear in q vanishes. Terms quadratic, cubic, etc. in q represent changes 
in frequency, anharmonicity, etc. of that bath mode following optical excitation of chro- 
mophore 1. Because a collective coordinate of the type q will generally be of low frequency 
(long wavelength), its change in frequency, anharmonicity, etc. will be assumed to be negli- 
gibly small. The leading term in Eq. (^) is thereby taken to be sensibly g-independent as 
well as independent of the intramolecular coordinates; it can be absorbed into the nuclear 
Hamiltonian Heg. Hence, we arrive at a working model in which all optically active col- 
lective coordinates are treated explicitly, the remaining (optically inactive) coordinates are 
relegated to the bath, and a phenomenological treatment of electronic dephasing need not 
be introduced. The remaining g-dependence in Eq. ( |64|) resides in terms linear in (Qi — A^^) 
and {Q2 — Ag^). These g-dependent terms govern bath-induced vibrational relaxation in the 
site state \eg). 

In order to avoid unnecessary complications, we make an additional assumption that 
Qi and Q2 couple to orthogonal combinations of bath coordinates, denoted by gi and g2, 
respectively. The system-bath coupling becomes 

Veg{Qi^Q2.quq2) = V^l\Q^,q,) + V^^\Q2^q2) 

- ^^' " ^' ^ dQ[ + ^^' " ^' > dQ2 ■ ^ ^ 

There is an analogous expression for Vge{Qi,Q2,qi,q2)- Since we are neglecting frequency 
shifts of the bath modes, it is consistent to make linear approximations to the derivatives in 
Eq. di: 

oQi oqidQi 
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and similarly for — ^^S^tt — ^ • ^^ the end, assuming the bilinear coupling constants to be 

the same for both pairs of intramolecular and bath modes in both electronic excited states 
allows us to write simply 

H,,, = \eg){eg\{{Q, - Af )/,gi + {Q2 - Af)fiq,)} 

+ l^e)(^e|{(gi - Af )/,gi + {Q2 - AD/.g^)} • (68) 

In practice, we hold it that electronically exciting one chromophore leaves an intramolecular 
mode on the other undisplaced, so that Ag^ = Af^ = 0. 

5C. Vibronic Relaxation 

In Eq. ( pBD for the detector-frequency-weighted populations in the states \a')\b), the 
effects of the bath have been confined to the terms Tr6{(/?|e~*^o'^|7)(7'| pb e*^o'^|/?')}, which 
is the /5/3'*^ element of the reduced chromophore pair operator cr(t) obtained by propagating 
cr(0) = |7)(7'|Pfe under Hq and tracing over the bath. It is convenient to propagate the 
elements of the initial density matrix of the system term by term in order to facilitate the 
orientational averaging. Since the system-bath coupling is assumed to be weak, we can 
follow the recent work of Jean and co-workers and use Redfield theory to calculate the time 
development of the operators cr{t). 

Since Redfield theory is most naturally implemented here in the exciton (i.e. eigen- 
state) representation of the chromophore pair Hamiltonian, the system-bath coupling must 
be transformed from the site representation, in which the matrix elements of H^cp, Eq. (pBD, 
between site vibrational quantum number states are readily calculated, to the exciton repre- 
sentation. The system-bath coupling in the exciton representation connects various exciton 
states and ultimately leads to population and phase relaxation among the exciton states. 



The Redfield equations of motion for the reduced operators cr(t) have the form ||26| , |27|j29| 

0'f3t3'it) = -i£l3f3'Crf3f3'{t) + ^ Rpp' ,nK' (^nK'{t) (69) 

and the Redfield tensor Rpfj^nn' is specified by 
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where [H,,^{t)]p^ = e'''^'[H,,p]p^e-'^^' and (■■■)= Tn{- --p,}. 

The bath-chromophore pair interaction in Eq. (|68|) is a sum of products of chromophore 
pair and bath operators, and the transformation of the system-bath couphng operator to the 
exciton representation leaves the bath operators untouched. Therefore, bath-chromophore 
pair interactions in the exciton representation can readily be expressed as a sum of products 
of chromophore pair and bath operators, and we adopt Pollard and Friesner's notation for 
the Redfield tensor which is appropriate to such a form p7 



2 

Rp[5',KK' = '^[{Ga)t,'f3'{Ga)f3n{{Qaa)pK + (©aa)^'/?'} 
a=l 

-5.>p> E(Ga)/3c(Ga)«.(el)5. - V E(G'a)«'c(G.)5^.(0a~a).'d , (71) 

« i 

where (see Eq. (|68|)) 

Gi = \eg){eg\{Qr - Af )/, + \ge){ge\Q,fi 
G2 = \eg){eg\Q2fi + \ge){ge\{Q2 - Af)/, , (72) 

and iQta)i3S are given by 

iet)ps = / dTe-^^^^-{qai±T)qa) . (73) 
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The expression for the Redfield tensor given by Pollard and Friesner |2^ also includes con- 
tributions from cross- correlations of different sets of bath coordinates, but since we assume 
that gi and q2 are completely uncorrelated, the terms involving cross-correlations between 
gi and q2 vanish. 

Neglecting the imaginary part of {Qta)i3Sy Eq. (^) can be rewritten as 

iQaa)0S = ^ / rfre-''^-(g.(±r)g„) (74) 

^ J— CO 

1 f^ 

= e^^'^'^- / dTe-''^'^{qaqa{±r)) . (75) 

^ J — oo 
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Equations (^) and (|75D can be combined to obtain {Qaa)i3S in terms of the symmetrized 
correlation function, 

1 roo 1 

(®-)/^^ = TT^±(^ J-oo dre-''''^-{{qaqai±r)) + (ga(±r)g,)} , (76) 

Pollard and Friesner have emphasized that whatever its shortcomings, the neglect of the 
imaginary part of {Qta)i3S does not undermine detailed balance. In fact, the symmetrized 
correlation function in Eq. (|76D can be replaced by any approximate function that is real and 
even with respect to time. We introduce a mean-squared fluctuation and bath correlation 
time by hypothesizing that 

li{qaqa{±r)) + {q^{±r)q^) = (g,^)e-l^l/- , (77) 

Using Eq. (^ in Eq. (|76|) leads to 

When Tc is short compared to all relevant bohr periods in the system, we obtain 

This panoply of assumptions and approximations allows us to specify the entire Redfield 
tensor by choosing a single parameter. ||2^ In the absence of excitation transfer, the Redfield 



tensor element governing the rate of population decay from \eg)\lj^g) (= |1)) to \eg)\Qe^g) 
(= |0)) is 

i?oo,ii = (Gi)io(Gi)oi{^^^ • 2tMI)} , (80) 



where (6*1)10 = (6*1)01 = fi/ \/^, and we can therefore choose the product of ff and Tc,{q\) 
that yields a specified vibrational population decay rate from \eg)\lj^g) to |ef7)|0e0c,). 

In our calculations we make the secular approximation, namely, we keep only those 
elements of the Redfield tensor Rppi^^K' in which the element is not significantly smaller than 
the frequency mismatch \ujj3p1 — tu^K'!- 
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5D. Anisotropy Decays 

Figure 7 shows calculated anisotropies for a chromophore pair with molecular parameters, 
laser pulse duration, and detector time and frequency resolution all the same as those used 
for the calculations of Figure 1. The initial state of the chromophore pair is \gg)\OgOg) and 
the bath is at 300 K. The population relaxation rate (in the absence of energy transfer) from 
\eg)\leOg) to \eg)\OeOg) is chosen to be 1.5 ps^^ and for times less than 150 fs (overlapping 
laser pulse and detection window) the dissipative effects of the bath are not included. This 
is done in order to properly include the contributions to the fluorescence rate present when 
the excitation pulse and the detection window overlap. We combine the instantaneous 
fluorescence rate obtained from an isolated chromophore pair for times less than 150 fs 
with the rate including dissipation for times longer than 150 fs. The combined rate is 
then integrated and used to calculate the time-resolved signal. Since the combined rate is 
composed of two different calculations, there is a discontinuity in the combined rate, leading 
to a slight discontinuity in the slope of the short-time anisotropy. 

The dissipative effects of the bath become most evident after about 0.75 ps, where 
the calculated anisotropies in Figures 1 and 7 cease to resemble each other. After 1.5 
ps, the absence of oscillations in the time-resolved anisotropies of Figure 7 suggests that 
the coherences between the exciton states are almost completely destroyed by the bath, 
and the slow decay of the anisotropies beyond 1.5 ps results from the populations in the 
reduced density matrix of the chromophore pair approaching their equilibrium values at 300 
K. Our claim that the slow monotonic decay of the time-resolved anisotropies beyond 1.5 
ps is due to the thermal equilibration of populations is supported by the calculations of 
the site excitation probability presented below. The lower (thin) curve in Figure 7 shows 
a more significant decay in the 1.5-4 ps region than the upper (thick) curve because the 
excitation conditions that result in the anisotropy shown by the lower (thin) curve prepare 
a higher-lying superposition of states that subsequently decays into states localized on (the 
lower-energy) chromophore 1. Excitation conditions resulting in the anisotropy shown by 
the upper (thick) curve prepare a superposition of states fairly well localized on the lower- 
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energy chromophore, and therefore, population relaxation does not result in a significant 
decay of the anisotropy in the 2-4 ps region. 

Figure 8 presents a calculation of the time-resolved anisotropy designed to check whether 
the time-resolved anisotropy behaves in a qualitatively similar manner when higher-lying 
thermally populated initial ground states are included. The contributions from the lowest 
three thermally populated levels \gg)\OgOg), \gg)\lgOg), and \gg)\Oglg) have been included in 
Figure 8. The contributions to the anisotropy from \gg)\lgOg) and \gg)\Oglg) are qualitatively 
similar to the \gg)\OgOg), and calculated anisotropics in Figure 8 show a similar behavior to 
those shown in Figure 7. In the case of excitation 100 cm^^ above the zero-zero transition 
frequency of the higher-energy chromophore, there is a perceptible smoothing of the coherent 
oscillations at early times, compared to its counterpart in Figure 7. 

The initial anisotropy in Fig. 8 still remains above 0.4 even with the inclusion of higher- 
lying thermally populated initial ground states. For example, consider the initial state 
\gg)\lgOg). The short-time excitation is dominated by the strongest Franck-Condon transi- 
tions and occurs to the states \ge)\lgOe) and \eg)\leOg). (Recall that for very short delays 
both the excitation pulse and the detection window are effectively shortened, increasing 
their effective bandwidths.) Because each of these singly-excited states involves excitation 
on a different chromophore, the emission from these states back to \gg)\lgOg) is the type 
that leads to a larger than 0.4 initial anisotropy. 

Figure 9 shows a calculation of the time-resolved anisotropy, including vibrational relax- 
ation, for the same case of a pair of identical chromophores and the same excitation and 
detection parameters used in Figure 4. The vibrational population decay rate of 1.5 ps^^ 
from \eg)\leOg) to \eg)\OeOg) in the absence of energy transfer was again used in this calcu- 
lation. The overall qualitative behavior of the two anisotropy curves is in some ways similar 
to their counterparts in Figure 7. However, the slow decay in the 2-4 ps range seen in the 
anisotropy curves in Figure 7 is absent in the corresponding curves in Figure 9. This is 
consistent with the fact that there can be no net population transfer between two identical 
chromophores. 
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The net site excitation probability difference Pi(t) — P2(t), given by Eq. ( p3D for the case 
of the isolated chromophore pair, can also be calculated with the inclusion of dissipative 
effects of the bath. The net site population difference can be written in terms of the reduced 
density matrix of the chromophore pair 

Pi{t) - P2{t) = Tr,p{{P, - P2)Tr4e-^^otp^p(0)p,e^^o*}} , (81) 

where H'q = Hep + Hh + Hbcp, and Trcp{- ■ ■} denotes the trace over the states of the chro- 
mophore pair. The initial density matrix of the chromophore pair, neglecting the effects of 
the bath during the excitation pulse, is given by 



Pcp(O) = J2 \l){l'\l^a'yl^ay{eaj ' eL)(e„y ■ eL)F{e^a - ^)F{eYc. - fi) (82) 

Calculation of -Pi(t) — -PjI^) for the system parameters and excitation conditions identical 
to those used for obtaining the anisotropics in Figure 7 are shown in Figure 10. The upper 
(thick) curve in Figure 10 is initially similar to the thick curve in Figure 6a, where the effects 
of the bath were not included. Pi(t) — P2{t) for the isolated chromophore pair is oscillatory 
but remains close to its initial value. When the bath is included, the oscillations are damped 
out and the excitation probability difference tends to its equilibrium value. The behavior 
of the excitation probability difference shown by the lower (thin) curve of Figure 10 is even 
more dramatic. The large oscillations present in its analog in Figure 6b are completely 
damped out by 2 ps, indicating that coherent net population transfer has stopped. Both 
curves in Figure 10 approach the same steady-state value of 0.34, consistent with thermal 
equilibration at 300 K. 

6. CONCLUSIONS 

We have presented calculations of time-resolved fluorescence anisotropy from a chro- 
mophore pair undergoing excitation transfer following excitation with an ultrashort laser 
pulse. An intramolecular mode for each chromophore was explicitly included in our model, 
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in order to understand the effects of vibrational structure on the anisotropy. We have in- 
cluded a model detector in order to take realistic account of the limitations imposed by the 
time and frequency resolution of the experimental detection process. 

Choosing parameters for our model system to correspond to the degree possible to known 
parameters in closely-coupled chromophore pairs in C-PC, we have investigated the excita- 
tion wavelength dependence of the time-resolved anisotropy, observing behavior qualitatively 
similar to the experimental findings. 0] Also, using a simplified version of our model that 
can be treated analytically, we have shown that the initial values of the anisotropy are highly 
dependent on the ground state(s) into which the chromophore is emitting. 

The inclusion of vibrational structure is shown to play a major role in determining the 
form of the time-resolved anisotropy. Several effects are observed that are not present if each 
chromophore is treated as an electronic 2-level system without nuclear degrees of freedom. 
First, the presence of multiple vibrational levels in the excited state of each chromophore 
introduces excitation-wavelength dependence that is significantly different from that found 
in a pair of coupled two-level systems. Second, in a chromophore pair lacking vibrational 
structure, coherent energy transfer requires a pulse bandwidth in excess of the bare exciton 
splitting, which is greater than the site energy difference. On the other hand, when each 
chromophore has a ladder of excited vibronic states, preparation of a superposition of close- 
lying vibronic states with amplitude on both chromophores can sometimes be accomplished 
with a pulse bandwidth less than the difference between the two sites in zero-zero transition 
frequency. Finally, the vibrational levels in the ground electronic state selected by the 
detection bandwidth play a significant role in determining the effects of quantum mechanical 
interference on the time-resolved anisotropy. 

We have reported Redfield theory calculations of the anisotropy, which included cou- 
pling of the intramolecular coordinates to a surrounding medium. When the effects of a 
thermal bath were included, the coherent oscillations in the anisotropy remained at short 
times but were seen to decay prior to the anisotropy reaching its steady-state value. It is 
interesting to note, however, that evidence for the coherent oscillations, which are present 
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in our calculations including vibronic relaxation, seems to be lacking in the experimental 
anisotropics measured for APC and C-PC by Xie et al. |^ Failure to sum over all vibra- 
tional levels initially thermally populated at 300 K causes the coherent oscillations in the 
anisotropy to be exaggerated. It can also be anticipated that addition of multiple optically 
active intramolecular and/or collective modes would tend to obscure the oscillatory contri- 
bution to the calculated signals. Furthermore, the coherent oscillations in the anisotropy can 
be degraded by the presence of inhomogeneous broadening in the zero-zero frequency offset 
between the two chromophores. Even if the zero-zero transition energies in the two chro- 
mophores are perfectly correlated, the inhomogeneous broadening of the zero-zero transition 
energies will have the effect of creating a distribution of different excitation and detection 
conditions, which could also tend to mask the coherent anisotropy oscillations. 

In very recent work, Jean [^ has reported calculations on fluorescence emission from an 
electron transfer system including vibronic relaxation via Redfield theory and a treatment 
of the frequency-dependent instantaneous emission rate essentially equivalent to our model 
detection apparatus. He reports the interesting result that under some circumstances, the 
instantaneous emission rate can accurately reflect the extent of electron transfer. 
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APPENDIX: ORIENTATIONAL AVERAGING 

We need to obtain the average over the products of four direction cosines in Eq. (^TJ), 
where e_|_ and e_ are unit vectors in the molecule-fixed frame, e^ and e^ are unit vectors in 
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the laboratory-fixed frame, and tliere is a constant relative angle 7 maintained between e^ 
and e_. The averages of the first two products in Eq. (|3T|) are |^ 



(e± ■ eL)^(e± ■ 6^)2 = - for bl \\ e^ 

= — for Bl Lea. (83) 

The average over the crossterm (e_|_ ■ ei)(e_ ■ e/^)(e+ ■ e(i)(e_ ■ e^^) can be calculated by 
expressing e_ in terms of e+ and e^^, a vector perpendicular to e^. 

e_ = e+ cos 7 + e_(._L sin 7 . (84) 

Substituting Eq. (|^) for e_ enables us to express the above product as a sum of products 
where the molecule-fixed unit vectors are either parallel or perpendicular to each other. The 



average can then be calculated to give [31 



(e+ ■ e^) (e_ ■ e^) (e+ ■ e^) (e_ ■ e^) = - cos^ ^ + 77 ^^^^ ^ ^^^ e^ || e^ 

= — cos^ ^ ~ ^n ^^^^ ^ f"^ eL -Lea . (85) 

Eq. ( pSD is obtained by noting that for the case of identical chromophores, we have 
/zi = /i2 = /i and e+ ± e_, (i.e. 7 = vr/2). Therefore, using Eq. (p6|), we can write 

fil = fj.^0g0g\0e0g)^il±ei-e2). (86) 

Defining cos0 = 61-62, and substituting Eq. ( ^61) into Eq. (pTf), Eq. (|38| ) is obtained. 

The orientational averages required for the evaluation of Eqs. (^9D and (^) can be 
similarly evaluated by writing 

62 = ei cos (j) + ei_L sin , (87) 

which makes it possible to express all orientational averages in terms of products of aver- 
ages involving only parallel or perpendicular unit vectors in the molecule-fixed frame. The 
averages needed in Eqs. ( ^Q]) and (p3| ) are 
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(ei(2) ■ 6^)2(61(2) • erf)2 = - for e^ || e^ 



15 



for bl ± erf 



1 .2^ 1 ... 



(62(1) ■ eL)^(ei(2) ■ 6^)2 = - COS + — sin for Bl \\ e^ 

1 2 

= — cos^ + 77 sin^ /or ez, ± e^ , 
15 15 



1 



1 



(62(1) ■ 6L)(ei(2) ■ 6L)(ei(2) " 6^) (62(1) ■ e^) = ^ cos^ *^ "^ 15 ^^'^'^ ^ ^^^ ^^ " ^'^ 

— cos^ 0-777: sin^ /or e^ ± e^ , 
15 3U 



and 



(62(1) ■ 6L)2(ei(2) ■ 6rf)(e2(i) ■ Gd) = - COS f Or G^ || Grf 



15 



COS /or gl ± erf 



(88) 
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t (ps) 

Figure 1. Time- resolved anisotropy plots for two different excitation/detection condi- 
tions. Upper (thick) curve corresponds to excitation at zero-zero transition frequency of 
lower-energy chromopliore, lower (thin) curve is for excitation 100 cm^^ above the zero- zero 
transition frequency of higher-energy chromophore. Detector, with a time window of 70 fs, 
is red-shifted from excitation by 250 cm^^. 




t (ps) 

Figure 2. Time-resolved anisotropy corresponding to excitation at 100 cm^^ above zero- 
zero transition frequency of higher-energy chromophore. Thin curve results from time de- 
tection window of 70 fs (same as in Fig. 1), and thick curve is for a time detection window 
of 150 fs. 
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t (ps) 

Figure 3. Time-resolved anisotropy for chromophore pair with 100 cm^^ difference be- 
tween the zero-zero electronic transition frequencies. Upper (thick) curve corresponds to 
excitation at zero- zero transition frequency of lower-energy chromophore, lower (thin) curve 
is for excitation 150 cm~^ above the zero-zero transition frequency of higher-energy chro- 
mophore. 




t (ps) 

Figure 4. Time-resolved anisotropy for a pair of identical chromophores. Upper (thick) 
curve corresponds to excitation at zero- zero transition frequency of the chromophores, lower 
(thin) curve corresponds to excitation 250 cm^^ above the zero-zero transition frequency. 
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Figure 5. Anisotropy at t = as a function of excitation transfer coupling. Excita- 
tion/detection conditions and model system parameters are the same as in the lower (thin) 
curve of Figure 1. 




Figure 6a. Thick curve shows the site population difference obtained with same system 
parameters and excitation conditions as in the upper (thick) curve of Fig. 1. The thin curve 
shows the corresponding anisotropy. Both ^ },{ n ),{ and R(t) are dimensionless. 
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Figure 6b. Thick curve shows the site population difference obtained with same system 
parameters and excitation conditions as in the lower (thin) curve of Fig. 1. The correspond- 
ing anisotropy is shown with a thin curve. 




Figure 7. Time-resolved anisotropy for chromophore pair with same parameters as in 
Figure 1 and a 1.5 ps^^ vibrational relaxation rate. Upper (thick) curve is for excitation at 
zero-zero transition frequency of lower-energy chromophore, lower (thin) curve is for excita- 
tion 100 cm~^ above the zero-zero transition frequency of the higher energy chromophore. 
The vibrational and electronic ground state of the chromophore pair is the only initial state 
included. 
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Figure 8. Time- resolved anisotropy for chromophore pair with the same parameters as 
in Figure 1 and a 1.5 ps~^ vibrational relaxation rate. Three lowest thermally populated 
initial states have been included. 



R(t) 




Figure 9. Time-resolved anisotropy for pair of identical chromophores with the same 
parameters as in Figure 4 and a 1.5 ps^^ vibrational relaxation rate. Upper (thick) curve 
corresponds to excitation at zero-zero transition frequency of either chromophore, lower 
(thin) curve is for excitation 250 cm~^ above the zero-zero transition frequency of either 
chromophore. 
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Figure 10. Site population difference for the chromophore pair with 150 cm^^ offset 
between the zero-zero electronic transition frequencies and a 1.5 ps^^ vibrational relaxation 
rate. Upper (thick) curve is for excitation at zero- zero transition frequency of lower-energy 
chromophore, lower (thin) curve corresponds to excitation 100 cm^^ above the zero- zero 
transition frequency of higher-energy chromophore. 
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